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Abstract 

The purpose of this work is to use a renormalized quantum scalar field to investigate very early 
cosmology, in the Planck era immediately following the big bang. Renormalization effects make 
the field potential dependent on length scale, and are important during the big bang era. We use 
the asymptotically free Halpern-Huang scalar field, which is derived from renormalization-group 
analysis, and solve Einstein's equation with Robertson- Walker metric as an initial-value problem. 
The main prediction is that the Hubble parameter follows a power law: H = a/ a ~ t~ p , and the 
universe expands at an accelerated rate: a ~ exp t l ~ p . This gives "dark energy" , with an equivalent 
cosmological constant that decays in time like t~ 2p , which avoids the "fine-tuning" problem. The 
power law predicts a simple relation for the galactic redshift. Comparison with data leads to the 
speculation that the universe experienced a crossover transition, which was completed about 7 
billion years ago. 
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I. INTRODUCTION AND SUMMARY 



According to quantum field theory, the vacuum is not empty and static, but filled with 
fluctuating quantum fields. Those of the electromagnetic field, which fluctuate about zero, 
can be measured experimentally through the Lamb shift in the hydrogen spectrum, and 
the electron's anomalous magnetic moment. Others, such as the scalar Higgs field of the 
standard model, fluctuate about a nonzero vacuum field. Grand unified models call for 
still more vacuum scalar fields. These vacuum scalar fields are similar to the Ginsburg- 
Landau order parameter in superconductivity, which is a phenomenological way to describe 
the condensate of Cooper pairs of the more fundamental BCS theory. Be they elementary or 
phenomenological, these vacuum fields behave like classical fields in many respects. Under 
certain conditions, however, one must take into account their quantum nature. In particular, 
during the big bang, when the length scale of the universe undergoes rapid change, one 
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must take into account the effects of renormalization, and this is the focus of the present 
investigation. Some of our results have been reported in a previous note [1]. 

Scalar fields have been used in traditional cosmological theories to explain "dark en- 
ergy" [2], and "cosmic inflation" [3]. Dark energy refers to an accelerating expansion of the 
universe, which can be reproduced by introducing a "cosmological constant" in Einstein's 
equation. This is equivalent to introducing a static scalar field with constant energy density. 
The problem is that the cosmological constant is naturally measured on the Planck scale, 
which is some 60 orders of magnitude greater than that fitted to presently observed data. 
One would have to "fine-tune" it (by 60 orders of magnitude!), and this has been deemed 
unpalatable. 

The theory of cosmic inflation, designed to explain the presently observed large-scale 
uniformity of the universe, postulates that matter was created while the universe was so 
small that all matter "saw" each other. The universe then expanded by an enormous order 
of magnitude (e.g., 27) in an extremely short time (e.g., 10~ 26 s), pushing part of the matter 
beyond the event horizon of other parts, but the original density was retained. To implement 
this scenario, one introduces a scalar field with spontaneous symmetry breaking, i.e., having 
a potential with a minimum located at a nonzero value of the field. Initially the universe 
was placed at the "false vacuum" of zero field, and it is supposed to inflate during the 
time it takes to "roll down" the potential towards the true vacuum. It would be desirable 
to formulate this scenario in terms of a mathematically consistent initial-value problem. 
However, this has not been done so far. As we shall see, the universe does inflate in our 
model, but the "rolling" was anything but slow. 

Most previous works on vacuum scalar fields treat them classically, i.e., with fixed given 
potentials. In quantum field theory, however, the potential is subject to renormalization, and 
changes with the energy scale. This arises from the fact that there exist virtual processes with 
momenta extending all the way to infinity. The high end of the spectrum causes divergences 
in the theory, and in any case does not correspond to the true physics. To make the theory 
mathematically defined, the spectrum must be cut off at some momentum A, and this cutoff 
is the only scale parameter in a self-contained field theory. When A changes, all coupling 
constants must change in such a manner as to preserve the theory (i.e., to preserve all the 
correlation functions), and this change is called "renormalization". Such cutoff dependence 
can be ignored when one studies phenomena at a fixed length scale, such as stellar structure 
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at a particular epoch of the universe; but it is all-important during the big bang. 

The purpose of this work is to study the implications of renormalized quantum scalar 
fields in the immediate neighborhood of the big bang. The mathematical problem is to 
formulate and solve an initial-value problem based on Einstein's equation, with suitable 
idealizations to render the problem tractable. This basic principle is that there is only one 
scale in the early cosmos, namely the "radius" a of the universe set by the metric tensor. 
Thus, we must identify a with inverse cutoff momentum A -1 . For consistency, the self- 
interaction potential of the scalar field should be "asymptotically free", i.e., vanish in the 
limit a — > 0. 

From renormalizat ion-group (RG) analysis, Halpern and Huang (HH) [4] have shown that 
asymptotic freedom determines the potential of the scalar field to be a Kummer function, 
a transcendental function that has exponential behavior for large fields, and this rules out 
all polynomial potentials, including the popular 4 theory. In the present work, we use the 
HH scalar field as the source of gravity, in Einstein's equation with Robertson- Walker (RW) 
metric. As mentioned earlier, our basic principle is that the inverse radius of the universe 
acts as the momentum cutoff of the scalar field theory, i.e., A = a -1 . This gives rise to 
a dynamical feedback: the expansion of the universe is driven by the scalar field, whose 
potential depends on the radius of the universe. 

The main prediction of the model is that the Hubble parameter H = a/a behaves like 
a power H ~ t~ p (0 < p < 1), for large times, after averaging over small rapid oscillations. 
The exponent p depends on model parameters and initial conditions. This indicates "dark 
energy" , for the universe expands with acceleration, according to a ~ exp t x ~ v . This behavior 
corresponds to an equivalent cosmological constant that decays with time like t~ 2p , and this 
avoids the usual fine-tuning problem. The origin of the power law can be traced to a 
constraint on initial values from the 00 component of Einstein's equation. 

Although our model is valid only in a neighborhood of the big bang, it is hard to resist 
to compare it with observations from a much later universe. A partial justification for 
doing this is that the power-law character may survive generalizations of the model. In this 
spirit, we calculate the relation between luminosity distance di and red shift z for a light 
source, according to the power law. To an extremely good approximation, we find <1l{z) = 
z(l + z)d , in which the exponent p enters only through the constant do. Comparison 
with data on the galactic redshift, from supernova and gamma-ray burst measurements, 
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suggest that there was an epoch in which do had a different value from the current one, and 
connecting the two epochs was a crossover transition completed about 7 billion years ago. 

Finally we address the scenario of cosmic inflation, which is inseparable with matter cre- 
ation. The question is whether enough matter can be created for subsequent nucleosynthesis, 
during the time when the universe was small enough that all constituents remained within 
each other's event horizon. 

An equally important question has to do with the emergence of the matter energy scale, 
which is from the Planck scale by some 18 orders of magnitude. Physically, the matter scale 
arises spontaneously, via "dimensional transmutation" in QCD, and in our model it enters 
through the coupling parameter between matter and the scalar field. These two scales must 
decouple from each other. How does it happen mathematically in our model? 

To explore these questions, we treat matter a perfect fluid coupled to the scalar field, 
as detailed in Appendix C. Our studies lead to the opinion that a completely spatially 
homogeneous scalar field, real or complex, cannot give a satisfactory inflation scenario. 
First, it cannot created enough matter in a short enough time, and secondly decoupling 
does not occur, whatever one chooses for the matter coupling parameters. The model so far 
appears to lack physical mechanisms for matter creation and decoupling. 

We are led to investigate a complex scalar field with uniform modulus but spatially 
varying phase. This makes the universe a superfluid, and new physics emerges, namely 
vorticity and quantum turbulence. We find that these phenomena can supply the missing 
mechanisms for matter creation and decoupling. This development is the subject of paper 
II of this series [5]. 



II. PRELIMINARIES 

We start with Einstein's equation 

Rw - t;9^ r = SnGT^ (1) 

where is the metric tensor that reduces to the diagonal form (—1, 1, 1, 1) in flat space- 
time, is the energy-momentum tensor of non-gravitational fields, and G = 6.672 x 1CT 11 
m 3 kg -1 s~ 2 is the gravitational constant. We shall put AnG = 1, thus measuring everything 
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in Planck units [6]: 

Planck length = (he' 3 ) 1 ' 2 (4ttG) 1/2 = 5.73 x 1(T 35 m 
Planck time = (ftc" 5 ) 1/2 (AnG) 1/2 = 1.91 x 1(T 43 s 
Planck energy = (he 5 ) 1/2 (4nGy 1/2 = 3.44 x 10 18 GeV (2) 

Consider a spatially homogeneous universe defined by the Robertson- Walker (RW) metric, 
which is specified through the line element 

ds 2 = -dt 2 + a 2 (t) ( + r 2 d6 2 + r 2 sin 2 6d<A (3) 
\ 1 — kr 2 J 

where t is the time, {r, 9, 0} are dimensionless spherical coordinates, and a (t) is the 
length scale. The curvature parameter is k = 0,±1, where k = 1 corresponds to a space 
with positive curvature, k — — 1 that with negative curvature, and k = is the limiting case 
of zero curvature. With the RW metric, the 00 and ij component of Einstein's equation 
reduce to the following Friedman equations: 

d\ 2 k 2 



- H — o — ~ o-^ 00 
a I a z 3 



d f d\ 2 k 

-+(-} +~ 

a V a I a z 



g tJ = -2T lJ (i,j = 1,2,3) (4) 

It is customary to introduce the Hubble parameter defined by 

H= h - (5) 
a 

The energy-momentum tensor of a spatially uniform system must have the form 

T 00 =p 
T ij = g ij p 

T j0 = (6) 

where p defines the energy density, and p the pressure. Energy-momentum conservation is 
expressed by T?" = 0, which, with the RW metric, becomes 

P + ^( P + P)=0 (7) 
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We can recast the Friedman equations in terms of H, and, with inclusion of the conservation 
equation, obtain three cosmological equations: 

H = - 2 -(p + p) 

a z 

TT o k 2 
p=-3H(p + p) (8) 
The second equation is a constraint of the form 

* sfl . + *_|, = „ (9 ) 

The third equation, the conservation law, states X — 0, i.e., the constraint is a constant of 
the motion. 

As an example, consider Einstein's cosmological constant A , which appears in a static 
energy-momentum tensor of the form (with units restored for convenience) 

T^ u = -g^ (Ao/SttG) (10) 

Corresponding to this, the energy density and pressure are given by 

po = Ao/87rG 

Po = -A /8nG (11) 

The conservation equation now states po = 0, which is trivial. Thus the cosmological 
equations reduce to 

k 



H 

a 2 

k 2 



// 2 = -^ + §Po (12) 
The asymptotic solution describes an exponentially expanding universe, with 

a (t) ~ exp (Hoot) 

H^ = (A o /127rG0 1/2 (13) 

Since a(t) is accelerating, we can say that there is "dark energy". However, the "natural" 
value of Hoo should be of order unity on the Planck scale, whereas the presently observed 



Hubble parameter is of order 1CT 60 . One would have to "fine tune" H^, by sixty orders of 
magnitude. 

With a dynamical scalar field, the constraint implies = 0. This is illustrated in 
Appendix A in an exact solution for the massless free scalar field, in which a/ a decays 
according to a power law, which is equivalent to saying that decays like a power. The 
effective cosmological constant is being "fine-tuned to zero", so to speak. This "automatic 
fine-tuning" also happens in our model, to be discussed later. 

III. HALPERN-HUANG SCALAR FIELD 

The HH scalar field that we use in this work has an asymptotically free potential, which 
is summarized here. Appendix B give a derivation from renormalization theory. 

For generality, consider an A^-component real scalar field n (x) with O(N) symmetry, 
with Lagrangian density (with h = c = 1) 



where 2 = Yl n =i ^n- The high-energy cutoff A is introduced through a modification of the 
two-particle propagator at small distances. (See Appendix B for details.) The form of the 
modification is not important here; what is important is that A is the only intrinsic scale of 
the scalar field. All coupling constants g n in the power-series V = J2 n 9n4> n must scale with 
appropriate powers of A. In 4-dimensional space-time we have g n = A A ~ n u n , where the u n 
are dimensionless, but depend on A; they undergo "renormalization" in order to preserve 
the theory. As A changes, {u n } trace out an RG trajectory in parameter space. There exist 
fixed points in this space, representing scale-invariant systems with A = oo. A obvious fixed 
point is the Gaussian fixed point corresponding to V = 0, i.e., the massless free field. 
In a universe governed by the RW metric with length scale a, we must identify 

A = 5 (15) 

where we restore Planck's constant h to remind us of the quantum nature of the cutoff. The 
big bang corresponds to a = 0, or the Gaussian fixed point. In a consistent theory, therefore, 
the potential must vanish as a — > 0, or A — > oo. In the language of particle physics, the 
theory must be "asymptotically free". We imagine that at the big bang, the scalar field was 




(14) 



n=l 
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displaced infinitesimally from the Gaussian fixed point onto some RG trajectory, along some 
direction in the parameter space. This initial direction determines the form of V. If the 
trajectory corresponds to asymptotic freedom, i.e., if the Gaussian fixed point appears as an 
ultraviolet fixed point on the trajectory, the potential will grow to engender a universe. A 
trajectory that is non-free asymptotically is a critical line on which all points are equivalent 
to the fixed point, and the system behaves as if it had never left the fixed point, with the 
time development as described in Appendix A. 

All quantities with dimension scale with appropriate powers of A. The potential V is of 
dimensionality (length) -4 , and we introduce a dimensionless potential U by writing 

V = A 4 U (16) 

Under a scale transformation, U changes under renormalization according to 

dU 

= P [U] (17) 

where the "beta-function" /3[U] is a functional of U. Near the Gaussian fixed point, where 
U = 0, we can make a linear approximation 

/3 [U] « -bU (18) 

leading to an eigenvalue equation 

= -m (i 9 ) 

which defines the eigenpotential U^. In the linear approximation, the most general U is a 
linear superposition of these eigenpotentials. 

From the renormalizat ion-group analysis briefly summarized in Appendix B, one obtains 
the solution 



U b (z) = cA~ b [M (-2 + b/2, N/2, z) - 1] 

z = Stt 2 J2 (20) 

n 

where M is a Kummer function, c is an arbitrary constant, and {p n (x) is a dimensionless 
field: 

(x) = -^t) n (x) (21) 
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Again, we restore units to remind us that the potential depends on h. 

The power series and asymptotic behavior of the Kummer function are given by 



' q q(q+l)2\ q(q + l)(q + 2) 3! 

M(p, q, z) « T (g) T^ 1 (p) ^ exp z (22) 

Using the derivative formula [7] 

M'(p,q,z) =pq- 1 M(p + l,q+ 1, z) (23) 



we obtain 



I7£(z) = -cA^A^ 1 (4 - b) M (-1 + 6/2, 1 + iV/2, z) (24) 



Asymptotic freedom corresponds to 6 > , and spontaneous symmetry breaking occurs 
when b < 2 . Thus we limit ourselves to the range < b < 2. 

The limiting case 6 = 2 corresponds to the massive free field, which is asymptotically free 
but does not exhibit spontaneous symmetry breaking, i.e., it does not maintain a vacuum 
field. The limiting case 6 = corresponds to the 4 theory, which exhibits spontaneous 
symmetry breaking, but is not asymptotically free. In our linear approximation (fl9j) . 6 = 
corresponds to AdU/dA = 0, which indicates neutrality. However, the beta-function to 
second order gives [8] 

A M" = ilL^ 2 ( for b = °' or 04 theory ) (25) 

which shows it increases as A increases, and is thus asymptotically non-free. 

We emphasize that the HH eigenpotential is derived (a) in flat space-time, (b) in the 
neighborhood of the Gaussian fixed point, where {7 = 0. Corrections due to space-time 
curvature and nonlinearity in U have not been calculated; but the present approximation 
should be good in a neighborhood of the big bang. 
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IV. COSMOLOGICAL EQUATIONS 



The canonical Lagrangian (Tl4j) of the scalar field gives the following equation of motion 
and components of the energy-momentum tensor : 

dV 



(j) n = -3H<f) 

1 N 

Pcanon = ~ ^ <t>n + V 
n=l 

N 

2 



71=1 

The constraint equation now reads 



X = H2 + yiY,ti-lv = (27) 



On general principle, the equations of motion must guarantee X = 0, since it is known that 
the Cauchy problem in general relativity exists [9]. However, direct computation using X 
as given in f )271) yields X = — (2/3) a (dV/da), which is nonzero if the cutoff depends on the 
time. This defect can be attributed to the fact that the gravitational cutoff has not been 
built into the Lagrangian (fl4|) . As remedy, we modify T^ u of the scalar field by adding a 
term to the pressure, and take 

P Pcanon 

adV . . 

P = Pcanon - g (28) 

For an eigenpotential V = a~ A Ub it can shown that 

a _ H6 _4)V +I > (29) 

n 

The cosmological equations now become 



• k ^ - 2 1 dV 



3 da 



-3H(f> r , 



X = H 2 + 



dV 



-V = 



(30) 
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The first two equations now imply X = 0, and we have a closed set of self-consistent 
equations. 

We are able to work with a set of classical equations, because we have neglected quan- 
tum fluctuations about the vacuum scalar field. However, important quantum effects are 
incorporated through the scale dependence of the potential V arising from renormalization. 

V. CONSTRAINT EQUATION AND POWER LAW 

The constraint equation in (l3"0j) requires 




That H be real and finite imposes severe restrictions on initial values. In particular, a = is 
ruled out; the initial state cannot be exactly at the big bang. This poses no problem from a 
practical point of view, for an initial universe with radius a ~ 1 (Planck units) is practically 
a point. 

From a physical point of view, we do not expect the model to be valid in the immediate 
neighborhood of the big bang, which would be dominated by quantum fluctuations. The 
universe could have been created at very high temperatures, and rapidly cooled through a 
phase transition to reach a vacuum with spontaneous broken symmetry. Or it could have 
been be created in the broken state. There is no way to know what actually happened; all 
we know is that we start our model at some time after the big bang, but still in the Planck 
era, with a vacuum field already present. 

Now we turn to the consequence of the constraint. Since V = a~ 4 U, it would vanish 
rather rapidly in an expanding universe. The same is true of <f) n , which is proportional 
to a -1 by dimension analysis. Thus, the constraint ( l3Tj) would make H — > 0. Given the 
absence of relevant scale, we expect H to obey a power law: 

H ~r p 

a ~ exp t l ~ p (32) 

The argument is far from rigorous, of course, but the result is support by the exactly 
solution for the massless free field (Appendix A), and is verified in numerical solutions to 
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be discussed. The latter show that the power law emerges after averaging over small high- 
frequency oscillations. 



VI. NUMERICAL SOLUTIONS 



For numerical solutions, we limit ourselves to the simplest real scalar field (N = 

1). A multi-component field would yield qualitatively the same results for a completely 
uniform universe. It is convenient to rewrite the cosmological equations as a set of first- 
order autonomous equations: 



a = Ha 

k 2 1 dV 
cr 3 oa 

<j> — v 

*-- < 33 > 

There are 4 unknown functions of time: a,H,<j),v. The initial values must be real, and 
satisfy the constraint 



2„ 1 ; , k 



1/2 



H = \o V + a<f> ~ T (34) 



3 ' ' 3 r a 2 

Although this relation is preserve by the equations, numerical procedures tend to violate it, 
and it is difficult to extend time iterations indefinitely. As an exploratory investigation, we 
have not looked into algorithm improvement. 

For completeness, we restate the HH potential V, which is generally a linear superposition 
of eigenpotentials V b : 

V b (<f>) = a- 4 U b (z) 

U b (z) = ca b [M (-2 + 6/2, 1/2, z) - 1] 

z = 8vrV0 2 (35) 



where M is the Kummer function. Some useful formulas are 



13 



k 


b 


c 


«o 


00 


00 


#0 


P 


-1 


1 


0.1 


1.00 


0.01 


0.1 


1.00 


0.81 





1 


0.1 


1.85 


0.17 


0.2 


0.91 


0.65 


1 


1 


0.1 


1.85 


0.19 


0.2 


1.70 


0.15 



TABLE I. Computation data: k = curvature; b,c = potential parameters; others = initial data; p 
= output exponent. 



9Vb fh /AT/ I M 

a— = (b - A)V b + 



da <90 
16n 2 a- 2 (f)Ul 



dV b _ 1C 2 -2 



U' b {z) = -c (4 - b) a b M (-1 + 6/2, 3/2, z) (36) 
The model parameters are 

Curvature: k — 1, 0, —1 

Eigenvalue: < 6 < 2 (37) 

Potential strength: c 

A pair of values {b, c} should be specified for each eigenpotential in V. The c's should be 
real numbers of either sign, such that V be positive for large (f>, and have a lowest minimum 
at ^ 0. 

First we use an eigenpotential with 6 = 1, which is shown in Fig.l at a = 1. As the 
universe expands, it will increase uniformly by a factor a{t). This property is a linear 
approximation that holds for sufficiently small a(t). Fig.2 shows numerical results for this 
potential, for curvature parameter k = 0. We see that H (t) oscillates about an average 
behavior consilient with a power law H ~ t~ p , with p = 0.65. The main source of uncertainty 
in p arises from the limitation on time iterations, due to numerical violation of the constraint. 
Numerical results for p from a number of runs are tabulated in Table I. 

Next we consider a superposition of two eigenpotentials: 

U (z) = ciU bl (z) - c 2 U b2 (z) 

^ = 1.6, ci = 0.1 (38) 
b 2 = 0.4, c 2 = 5.0 
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The locations ±z min and the depth U m \ n of the minima are functions of a, and are plotted 
in Fig.3. Because of the large ratio c<ijc\ = 50, U m \-n suddenly jumps at a near-critical value 
a c w 5. For a < a c , the minima of the can be approximated by two symmetrically placed 5- 
functions; the scalar field becomes trapped at values ±0i corresponding to the minima, and 
the model approaches the Ising spin model. Results of numerical solutions are shown in Fig. 
4, with curvature parameter is k — 0, and the initials conditions are ao = 1, 0o = 0, <po = 0.1. 

Figs. 2 and 4 show that the scalar field oscillates during cosmic expansion, contrary to 
the " slow- roll" picture of inflation. Closer examination show that the oscillation amplitudes 
are so large as to sample the exponential region of the potential wall. That is, the distinctive 
part of the HH potential, which makes it asymptotically free, plays an important role in the 
expansion of the universe. 



VII. COMPARISON WITH OBSERVATIONS 



Our model is valid only in the Planck era, and does not contain matter apart from the 
vacuum scalar field. We shall nevertheless compare the model with present observations, 
assuming that the power law H (t) ~ h t~ p will persist in the real universe. The index p 
depends on model parameters, which might change with conditions in the universe such as 
the temperature. For our analysis, however, we take p to be a constant. All quantities are 
measured in Planck units, unless otherwise specified. 

The age of the universe t and the present value H now = H (t ) are taken to be 

to = 1-5 x 10 10 yrs « 10 60 
H now = t (39) 

The initial value, defined at t — 1, is given by 

# initial = ho (1.65 x 10 50 )^ P) (40) 

If we put -^initial = 1 as a natural value, then h gives the fine-tune factor, which are tabulated 
Table I. 

The radius of the universe expands according to 

a it) = ao exp 

1 -p 



15 



V 


h 


0.5 


1.25 x 10 25 


0.85 


3 x 10 7 


0.95 


300 


0.99 


3 



TABLE II. Fine-tune factor for Hubble's parameter 



p 


Onow/ ^0 


0.5 


7.4 


0.85 


786 


0.95 


5 x 10 8 


0.99 


3 x 10 43 



TABLE III. Present radius of universe 

The present radius is a(l) : 

a now = a exp — !— (41) 
1 -p 

Some values are tabulated in Table II. 

Under the assumption that p is constant, its most reasonable value would lie in the range 
0.99 <p < 1. 

We now turn to the galactic redshift. The relation between the luminosity distance di of 
the source and the redshift parameter z is implicitly given by the following relations [10]: 



/(n) 



a (to) _ 1 
a(ti) 
t0 dt 



ti 



a(t) 



d L = "-^jj^- = na (t ) (1 + z) (42) 

where t the the time of detection, at the origin of the coordinate system, of light emitted 
at time t\ < t , by a source located at co-moving coordinate r\. The function / is defined 
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by 

sin -1 T\ ik = 1) 



f n _dr_ 

/(ri)=/ Tra 



ri (k = 0) (43) 
sinh -1 r\ (k = — 1) 



Using the first two equations, we can expressed r\ and t\ in terms of to and z, and then 
obtain di (z) from the third equation. 

In our model, a(t) = a exp (£t 1_p ) , where £ = /i (1 — p)~ ■ Define an effective time 
r = For < p < 1, the second equation in ([12]) can be rewritten as 

/"TO 

/ ( ri ) = K / rfrr p/(1 - p) exp (-r) (44) 
where K = [(1 - p) a ] _1 ^-Va-p), an d 

ro = ^ 

n = r -ln(z + l) (45) 

Since t ~ 10 60 , we can assume r >> 1, and obtain to a good approximation / (r x ) £s 
-ft'i-z, where K\ = Kqt^ 1 p ^ exp (— to). Since Kq is extremely small, this gives r\ = z to a 
very good approximation, and thus 

d L = JTia « (1 + (46) 

We rewrite this as 

^ = d r/ (1 + z) (47) 

where d = c/H now = 4283 Mpc, corresponding to the choice H now = 70 km s^Mpc -1 . 

Fig. 5 shows comparison with data from observations and supernovas [11] and gamma-ray 
bursts [12]. The upper panel shows the parameter \x used in conventional data analysis: 

p = h\o%[-^\ + 25 (48) 



Mpc 

plotted as a function of z. The lower panel shows a semilog plot of d^/z vs. z. Lines 
corresponding to Hubble's law (no dark energy) are shown. The p-dependence affects only 
the vertical displacement but not the shape of the model curves. Curve A corresponds to 
([47]) with 7] = 1, and curve B with 77 = 1/4. Curve A fits the data for z < 1, while curve B 
could represent the situation in a large- z regime beyond present measurements. 
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The power-law model allows only for variations in do, which may come from variations 
in the exponent p, caused by conditions such as the temperature. This leads us to specu- 
late that the universe may have had gone through a broad phase transition, or crossover, 
connecting two situations corresponding respectively to the curves A and B. The transition 
was completed around z — 1. 

The relation between the emission time and the red shift can be obtained from (j45p : 

^ = [1 - (1 - p ) In ( z + 1)] 1/(1 - p) (49) 
For p fa 1, we put p = 1 — e and obtain 

^[l-eln^ + l)] 1 ^— (50) 

to e->o 

Assuming this relation, we judge that the transition was completed at ti/t fa 0.5, or about 
than 7 billion years ago. 



VIII. COSMIC INFLATION AND DECOUPLING 



The problem of cosmic inflation is inseparable from that of matter creation, which has 
not been taken into account in our model so far. Most of the matter in the universe should 
have been created by the end of the inflation era, in order that the memory of the original 
density be imprinted. 

An equally important problem relates to energy scales. Our equations so far contains only 
one scale, the Planck scale. With matter creation, there emerges the scale of nuclear inter- 
actions, which is smaller than the Planck scale by some 18 orders of magnitude. Physically 
speaking, the matter scale emerges in QCD spontaneously through "dimensional transmuta- 
tion" [13]. In our model, we would introduce it through the coupling parameter between the 
scalar field and matter. These two scales must eventually decouple from each other. That 
is, the cosmological equations should break up into two sets, one governing the expansion, 
the other galactic evolution, and in each set the information about the other set occurs only 
through lumped constants. What is the mechanism for this decoupling? 

To address these questions, we model matter as a perfect fluid coupled to the scalar field, 
and obtain a set of cosmological equations that, again, represent an initial- value problem. 
These are derived in Appendix C. Numerical studies of these equations, both for a real scalar 
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field and a complex scalar field, lead us to the conclusion that a completely uniform scalar 
field, even with more than one components, cannot create sufficient matter to satisfy the 
inflation scenario. More important, it cannot exhibit the decoupling desired. 

We are led to an attempt to relax complete uniformity, within the dictate of the RW 
metric. It is natural to consider a complex scalar field with uniform modulus, but spatially 
varying phase. The phase variation gives rise to superfluid velocity, with the attendant 
vortex dynamics. The universe then becomes a superfluid with vortex dynamics. New 
physics emerges, namely the growth and decay of a vortex tangle that fills the universe, 
signifying quantum turbulence. We find that this provides a framework for matter creation, 
and the decoupling of scales. 

In the extension of our model, the demise of quantum turbulence will signifies the end 
of the inflation era, as well as the validity of our model, for density fluctuations would 
become important. The standard hot big bang theory will then take over, with one addition: 
the universe remains a superfluid with vorticity. The latter will offer explanations to post- 
inflation phenomena such as galactic voids, galactic jets, and the dark mass. We will present 
this development in detail in paper II of this series [5]. 



Appendix A: The massless free field 

The cosmological equations with a real massless scalar field, corresponding to V = 0, are 



a = Ha 
k 



i,2 



H =~- 

a 1 

= -3#0 

X^ 2 -i0 2 + A = O (Al) 

O CL 

They describe what happens if the scalar field remains at the Gaussian fixed point. The last 
equation X = is the constraint equation, and X is a constant of the motion. 
The third equation can be rewritten in the form din (j)a 3 ^j /dt = 0, which gives 

= c a~ 3 (A2) 
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where Co is an arbitrary constant. The equations then reduce to 

a = Ha 

H = — - ^2- 

a 2 a 6 

2 ci fc 



# = "T - "a ( A3 ) 



where ci = Cg/3. Dividing the second equation by the first, and equating H/a = dH/da,we 
obtain 

HdH = (^--%] da (A4) 



Integrating both sides gives 



Since H = a/a, this can be further integrated to yield 

/(ia . , . 

/ -4. 2 ». ( A6 ) 
VCia 4 + c 2 a 2 — fc 

where c 2 is an arbitrary constant. The ± signs reflect the time- reversal invariance of the 
equations. We choose the positive sign to obtain 

a (t) — y a exp (y/c^t) (A7) 



t— ¥OC 



This is the general solution without constraint, and c 2 is the equivalent cosmo logical con- 
stant. 

The constraint equation can be put in the form 



- = ±Vcia" 6 - ka- 2 (A8) 
a 

which gives 

/da , L . 

1 _ 4 r (A9) 
y/Cia 4 — ft 

Comparison with (??) shows 

c 2 = (A10) 

Thus, 0A7p is incorrect; the constraint "fine-tunes" the cosmo logical constant to zero. The 
correct solution gives 

' = c^V' 3 (k = 0) 
aft) I — ► q 1/4 (k = 1) (All) 
— >t (k = -l) 

k t— >oo 
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which corresponds to a power-law 



H — > hot' 1 

t^oo 



(A12) 



Appendix B: Renormalization and the Halpern-Huang potential 

A distinctive feature of quantum field theory is that the field can propagate virtually. 
This is described by the propagator function, which for a free field has Fourier transform 
A (A; 2 ) = k~ 2 . The high- A;, or high-energy modes must be cut off, for otherwise the virtual 
processes lead to divergences, rendering the quantum theory meaningless. The cut off energy 
A is introduced by "regulating" the propagator: 

f(z) -+ (Bl) 

z— s-oo 

The detailed form of f(k 2 / A 2 ) is not important. What is important is that A is the only 
scale in the theory. The regulated propagator in configurational space will be denoted by 
K (x, A) . 

In the formulation of renormalization according to Wilson [14,15], interaction coupling 
parameters must change with A, in such a fashion as to preserve the theory. This is called 
"renormalization". For a given value of A, the parameters define an effective theory appro- 
priate to that energy scale. A reformulation of the Wilson scheme using functional methods 
has been given by Polchinski [16]. 

Interactions that go to zero in the short-distance limit (or infinite-energy limit) are said to 
be asymptotically free, an example of which is the gauge interaction in QCD. In the opposite 
non-free behavior, the interactions grow indefinitely with decreasing length scale, and would 
diverge in the limit. This is the behavior found in QED and the 4 scalar field, for which 
the short- distance limit can exist only if there is no interaction at all. For applications in 
cosmology, we want interactions that vanish at the big bang, the small-distance limit, which 
means asymptotically free interactions. 

The Halpern-Huang (HH) potential was originally derived [4] by summing one-loop Feyn- 
man graphs. Here we outline an improved derivation due to Periwal [17], which is based on 
Polchinski's functional method of renormalization. For simplicity consider a real scalar field 
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(N = 1). The action in d- dimensional Euclidean space-time can be written as 



S[0,A] = S o [0,A] + S'[0,A] (B2) 

where the first term corresponds to the free field, and the second term represents the inter- 
action. We have 

S [0, A] = I J d d xd d y (x) K- 1 (x -y,A)<f> (y) (B3) 

where K^ 1 (x — y, A) is the inverse of the propagator K{x — y, A), in an operator sense. It 
differs from the Laplacian operator significantly only in a neighborhood of \x — y\ = 0, of 
radius A -1 . The partition function with external source J, which generates all correlation 
functions of the theory, is given by 

Z[J,A] = N J D0 e - s ^ A ]-W) ( B4 ) 

where M is a normalization constant, which may depend on A, and (J, 0) = J d d xJ (x) (x). 

In Wilson's renormalization scheme, modes contributing to the integral in (1B4j) with 
momentum higher than A are "integrated out", but not discarded, in order to lower the 
effective cutoff. This leads to a change the form of S', but the system itself is unaltered. 
The interactions are then said to be "renormalized" . In a general sense, renormalization 
means changing the cutoff A with simultaneous change in the form of S', so as to leave Z 
invariant, i.e., 

^ = (B5) 
dA 

This constraint is solved by Polchinski's renormalization equation, which is a functional 
integro-differential equation for S' [0, A]. For J = 0, it reads 



dS' If,, 9K (x — y, A) 

dxdy 



dA 2 J * OA [5(j)(x)5(f)(y) 5<p (x) 50 (y) 



S 2 S' 5S' 5S' 



(B6) 



Assuming that there are no derivative couplings, we can write S' as the integral of a local 
potential: 



5" [0, A] = A J d d xU (ip (x) , A) 

ip (:r) = A l - d,2 <p (x) (B7) 
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where U is a dimensionless function, and tp is a dimensionless field. In the neighborhood of 
the Gaussian fixed point, where S' = 0, we can linearize (IB6j) by neglecting the last term, 
and obtain a linear differential equation for U (if, A): 

A M + \ u " + G " t) vU ' + ud = (B8) 

where a prime denote partial derivative with respect to ip, and k = A 3 dK (0, A) /dA. Now 
we seek eigenpotentials Ub A) with the property 

A^ = -bU b (B9) 

In the language of perturbative renormalization theory, the right side is the linear approx- 
imation to the /3-function. Substituting this into the previous equation, we obtain the 
differential equation 

^ d 2 1/* ~\ d , ,\ 
W-2 {d - 2) ^ [d - h) 



u h = o 



(BIO) 



Since this equation does not depend on A, the A-dependence of the potential is contained 
in a multiplicative factor. In view of (1B9I) . the factor is A~ b . 
For d 2, (1B10I) can be transformed into Kummer's equation: 



where 



d 2 , . d 
dz z dz 



q=l/2 
b-d 



P 



U h = 



P 



d-2 

(2k)- 1 (d-2) if 2 



(Bll) 



(B12) 



The solution is 

U b (z) = cA' b [M (p, q, z) - I] (B13) 

where c is an arbitrary constant, and M is the Kummer function. We have subtracted 1 
to make Ub (0) = 0. This is permissible, since it merely changes the normalization of the 
partition function. In f !20p . the value of k corresponds to a sharp cutoff. 

For d = 2, the solution to (IBlOj) is sinusoidal, and the theory reduces to the XY model, 
or equivalently the so-called sine-Gordon theory [18]. 
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Appendix C: Coupling to perfect fluid 



We discuss how the cosmological equations (!30|) may be generalized to include coupling 
to galactic matter modeled as a perfect fluid, whose energy-momentum tensor is given by 
[19] 

T£ = -g^Pm + (p m + Pm) WIT (CI) 

where p m is the energy density, and is a velocity field, with g flu U tl U u = 1. For a spatially 
uniform fluid, U° = 1, W = 0. We assume the equation of state 

Pm = eoPm (C2) 

where €q = 1/3 for radiation, and eo = for classical matter. The coupling to the scalar field 
is specified via an interaction Lagrangian density £ mt . We give some examples of possible 
interactions. 

The simplest interaction is a direct interaction with a real scalar field: Ant = — Ap m 0. 
Current- current interaction with a complex scalar field (N = 2) can be constructed as 
follows. Represent the scalar field in terms of 4> = 2 -1 / 2 (</>i + ifo) and its complex conjugate 
0*, or in terms of the phase representation = Fexp(ia). The conserved scalar current 
density in the absence of interaction is J^ c = (2?)" 1 (<p*d^(f) — <f)d^<f)*) = F 2 d fl a. The current 
density of a perfect fluid is J™ = p m U u . The current- current interaction corresponds to 

Ant = -Xg^r; J? = Xfhng^F 2 (<V) U v 

= —Xp m F 2 & (spatially uniform system) (C3) 

Returning to the general case, we can decompose the total energy-momentum tensor of 
scalar field and perfect fluid as follows: 

= T s f + T£ + TZ (C4) 

We assume 

TZ = -<T A* (C5) 
which leads to an interaction energy density p\ nt and pressure p mt : 

Pint — — Ant 

Pint = Ant (C6) 
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The equation of motion for the perfect fluid comes from the conservation law Tf^ = 0, which 
for a spatially uniform system reduces to 



p + 3H(p + p) = (C7) 

where 



P = Psc + Pm + Pint = ^^4 ) l + V + Pm + £int 

n 

P = Psc + Pm + Pint = ^ J ^2'i ) l- V + e 0Pm ~ Ant 



(C8) 



We can rewrite (1C7I) in a more useful form. First, multiply both sides of the <p n equation 
in flU]) by n : 

<9V , «9£ int • 



<Pn<Pn = -3H(j) n - T—(pn + 



E<9£int ; 



(C9) 



We write 



Thus 



Using this we get 



dV _^dV ■ dV ■ 



dV TT 
+ a—H 

oa 



Now, using ()C8p . we can rewrite (IC7[) as 



(CIO) 



(Gil) 



(C12) 



d 
dt 



-3H 



^^ + (l + e ) 



Using the equation before this, we finally obtain 

dp m 



dt 



nrr /i j \ V- d£ int ; ^int ^ 

3H (1 + e ) p m - ~^I~ ( * >n ~ a —H 



dt 



da 



(C13) 



(C14) 
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In summary, the cosmological equations are, with H = a/a, 



H = A - 4ttG 
a 2 



1 OV 
H — Oj--— 
3 da 



!»„ = -SH(p n - — — + — — 




p m = -3if (1 + e ) p m - ^ - - ifa— 

n 

E^ + ^ + pJ-^ (C15) 

n=l / 

The last equation is a constraint on initial conditions, and is preserved by the equations of 
motion. This defines a self-consistent initial-value problem. 

Analytical and numerical studies show that matter creation is inefficient, and that no 
decoupling occurs between expansion and matter dynamics. 
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FIG. 2. Results from solving the initial- value problem with the potential U\ of Fig.l. The Hubble 
parameter H follows a power law t~ p after averaging over small oscillations. The flat tail is spurious, 
arising from numerical instability. The scalar field 4> oscillates with large amplitudes, sampling 
the exponential region of the potential. The behavior is quite different from the "slow-roll" of 
conventional inflationary theories. 
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min 




1CT 1 10° 10 1 10 2 10 3 

a 

a 

FIG. 3. The superposition of two eigenpotentials with a ratio of 50 in relative strength produces 
a potential with two symmetrically placed minimum that approach delta functions in the limit 
a — > 0. The scalar field becomes trapped in these minima, and the field theory approaches a spin 
Ising model. Here, the location of the minima ±z m ; n and potential depth J7 m i n are plotted as 
functions of a. 




I t Huang Fig.4 

FIG. 4. Results from solving the initial- value problem with superposition of eigenpotentials de- 
picted in Fig.3. 



30 




FIG. 5. Comparison between model prediction of the galactic redshift with observational data. 
Upper panel shows a conventional plot. Lower panel show a log-log plot of dj_,jz vs z, where di 
is the luminosity distance and z is the redshift parameter. The two theoretical curves, labeled A 
and B, correspond to different values of the exponent p, which depends on parameters in the scalar 
potential, and initial conditions. See text for fuller explanation. 
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